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Abstract 

Using coherent phase states, parameterized phase state distributions for a 
single-mode radiation field are introduced and their integral relation to the 
phase-parameterized field-strength distributions is studied. The integral ker- 
nel is evaluated and the problem of direct sampling of the coherent phase 
state distributions using balanced homodyne detection is considered. Nu- 
merical simulations show that when the value of the smoothing parameter 
is not too small the coherent phase state distributions can be obtained with 
sufficiently well accuracy. With decreasing value of the smoothing parameter 
the determination of the coherent phase state distributions may be an effort, 
because both the numerical calculation of the sampling function and the 
measurement of the field-strength distributions are required to be performed 
with drastically increasing accuracy. 



1 Introduction 



The quantum-mechanical description of amplitude and phase quantities and 
their measurement has turned out to be troublesome and is still of matter 
of discussion (for a review, see Ref. [0]). Attempts to introduce amplitude 
and phase operators in quantum mechanics and particularly in the theory 
of the quantized radiation field are nearly as old as quantum theory. Since 
Dirac's introduction of amplitude and phase operators in 1927 @, a series 
of concepts have been developed. Roughly speaking, there are two routes 
to introduce amplitude and phase quantities and the associated operators. 
In the first, amplitude and phase variables are introduced by intrinsic phase 
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closely related to phase space functions fjl 
second, they are introduced by defining them more pragmatically from the 



output observed in phase-sensitive measurements |22|, |23 



Phase measurements are usually based on some kinds of interference ex- 
periments. To our knowledge the earliest measurements of phase were car- 
ried out by Gerhardt, Buchler, and Litfin 
neered by Walker and Caroll [^5| 



A powerful method pio- 
for phase measurements has been homo- 
dyne and heterodyne detection techniques. These techniques, which were also 
used by Noh, Fougeres, and Mandel [23] for defining measurement-assisted 
phase quantities, have offered qualitatively new possibilities. In particular, 
Smithey, Beck, Faridani, and Raymer |26j have demonstrated experimentally 
that balanced homodyne detection can be used to reconstruct the quantum 
state of a signal- field mode (see also Ref. P8f ). The method is frequently 
called optical homodyne tomography. The Wigner function of the signal 
mode is reconstructed from the measured set of phase-parameterized field- 
strength distributions |27j] using inverse Radon transformation techniques, 
and the density matrix of the signal field is then obtained by Fourier trans- 
forming the Wigner function. It should be noted that tomographic methods 
have also been applied to the reconstruction of quantum states of molecular 
vibrations by Dunn, Walmsley, and Mukamel Since the Wigner func- 
tion or the density matrix contains all knowable information on the quantum 
statistics of the signal mode, they can be used for inferring phase distributions 
that have been not directly available from experiments, as was demonstrated 
by Beck, Smithey, and Raymer |30[ (see also Ref. [|3Tf ). 

The tomographic method for reconstructing the Wigner function and 
the density matrix of the signal mode from the phase-parameterized field- 
strength distributions measured in balanced homodyne detection (for the 
detection scheme see Refs. [[32|, ^4|) is based on the fact that knowledge 
of the field-strength distributions for all phases within a 7r-interval is equiv- 
alent to knowledge of the quantum state f35|. In other words, the measured 
set of field-strength distributions themselves can be regarded as being rep- 
resentative of the quantum state of the signal mode and the question arises 
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of whether or not quantum-phase distributions can be sampled directly from 
the measured field-strength distributions, without making the detour via the 
density matrix. The problem is similar to those related to the determination 
of the density matrix itself. Whereas in the first experiments the density ma- 
trix was obtained via the reconstructed Wigner function, recently methods 
of direct sampling of the density matrix from the measured data have been 



suggested and successfully applied [gg, \$7\, |38|, [39|, [40], ^T], |2|, [43], |44|, |45|] . The 



feasibility of measuring quantities by direct sampling depends on the relation- 
ships between the quantities and the field-strength distributions recorded. In 
particular, these relationships should be stable against sampling errors, so 
that measurement of a sufficiently large set of data should enables one to 
obtain the quantities with sufficiently well precision. Clearly, different con- 
cepts of introduction of phase variables lead to different phase distributions 
and to different relations between them and the field-strength distributions. 
In the present paper we restrict attention to the class of quantum-phase 
distributions that can be defined on the basis of the coherent phase states 



[p~0| , |I3"| , [I4|| . Roughly speaking, these distributions are smoothed London 
phase state distributions which can be be parameterized with respect to the 
degree of smoothing. The connection of the coherent phase state (CPS) dis- 
tributions with the phase-parameterized field-strength distributions is con- 
sidered and the problem of direct sampling of the CPS distributions from the 
field-strength distributions is studied. Numerical simulations are carried out 
to illustrate the applicability of the method. 

In Sec. |2| the CPS distributions are introduced. Their integral relation 
to the phase-parameterized field-strength distributions is studied in Sec. |[ 
Results of a computer simulation of direct sampling of CPS distributions are 
reported in Sec. [|, and a summary and some concluding remarks are given 
in Sec.||. 

2 Coherent phase state distributions 

In the long history of the quantum phase problem the eigenstates of the 
"exponential-phase" operator 

oo 

iL = ^|n)(n+l| (1) 

n=0 

(|n), photon-number states) have been played a central role. In particular, 
the states 

oo 

\e*+) = ]T e m *\n) (2) 

n=0 

satisfying the eigenvalue equation 

E_ | e **) = e <* \ e ^) (3) 
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have been referred absolutely to as phase states. These states introduced by 
London || resolve the unity, but they are non-orthogonal and non-normalizable. 
As was shown by Levy-Leblond P!D[ , the eigenvalue equation 

E_\z)=z\z) (4) 

has also other solutions than those in Eq. (El), viz. 



°° 

k) = v 1 - N 2 E Z 1 R )' M < L ( 5 ) 



n=0 



These normalizable states, which form an over-complete basis, have been 
re-examined and called coherent phase states by Shapiro and Shepard [|Hf|. 
Their properties have also been studied by Vourdas |14| who showed that 



they are SU(1,1) Perelomov coherent states \z;k) corresponding to the k = 
I representation [fi6 |. In the limit \z\ — > 1 the states \z) become the non- 



normalizable London states in Eq. (g): 

|e^)=lim^=L=|0,e), (6) 



,0 y/l 



where the notations 



e) = \z), z = e-V* (e > 0) (7) 



have been introduced. 

The states |0, e) can be used to define e-parameterized phase state distri- 
butions - CPS distributions - of a radiation-field mode via the overlap of its 
quantum state with coherent phase states: 

p(0,e)=iV- 1 (e)(0,e|^|0,e), (8) 

where 

N{e)= [ d<j>(<l>,e\Q\<j>,e) (9) 



(g, density matrix of a radiation-field mode). In the limit e — > the distri- 
butions p(<p, e) become the familiar London phase state distribution p(<p, 0) 
= (e l ^\g\e 1 ^) . The CPS distributions p(4>,e) can be regarded as smoothed 
London phase state distributions, which can be directly sampled from the 
field-strength distributions measured in balanced homodyne detection, as we 
will see below. 

To illustrate the effect of the smoothing parameter e, in Figs. [1] - |3] 
CPS distributions and the London phase state distributions of coherent and 
squeezed states are plotted for various values of the coherent amplitude and 
the squeezing parameter. When the value of the smoothing parameter is de- 
creased the CPS distributions become closer and closer to the corresponding 
London phase state distribution, which is observed, at least in tendency, in 
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the limit e — > 0. From comparison of Figs. |2] and |3] we further see that with 
increasing (mean) photon number the value of the smoothing parameter is 
required to be decreased in order to observe CPS distributions that are com- 
parably close to the London phase state distribution. This can be explained 
as follows. In the expansion in photon-number states of the London phase 
states the weight of the states with large photon numbers is increased with 
increasing (mean) number of photons in the state under study. In the CPS 
distributions the overlap of these states with the state under study would be 
suppressed when the smoothing parameter would be unchanged. 



3 Relation to the field-strength distributions 

The phase-parameterized field-strength distributions measurable in perfect 
balanced homodyning, 

p(f, <p) = (J 7 , <p\ Q \F, <p), (10) 

were introduced by Schubert and Vogel |27| (for details see, e.g., p7B). In 
Eq. fllCf ), the \ J-",<p) are the eigenvectors of a field-strength operator F(<p) = 
Fa + F*al at chosen phase ip (F= \F\e~ lip , a) and a are the photon creation 
and destruction operators, respectively). It should be noted that non-perfect 
detection introduces additional noise, so that the measured distributions are, 
in general, smeared field-strength distributions. These are convolutions of the 
true field-strength distributions with Gaussians p4| : 



P(F, <P\s) = J dFpiJ 7 ', if) v{T-T'\ s), (11) 

where 

V {T-s) = ( 2 ^Hi^lT lex p(-^^) ( 12 ) 

and s= 1 — r/ -1 , rj being the detection efficiency (rj < 1). 

The relation between the field-strength distributions and the CPS dis- 
tributions may be obtained as follows. Using the expansion of the density 



operator g in s-parameterized displacement operators |^7| an d relating the 
averages of the s-parameterized displacement operators to the characteristic 
functions of the field-strength distributions yields (see, e.g., Refs. |42|j) 

^ = So dip I ^P^i W s ) V 3 ; s), (13) 
where the operator integral kernel K^, tp; —s) reads as 

K(F,tp;-s) = 12 | d y\y\ exp{n/ [P(<p)-J=] - \sy 2 \F\ 2 }. (14) 
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Note that Eq. ( |HD can formally be applied also to cases where the detection 
efficiency is less than unity. Combining Eqs. (§) and (0), we find that p(4>, e) 
can be related to p(F, <p) as 

/*7T /*00 

p(0, e) = N-\e) / dip / dFp{F, if; s) K e (cf>, T, if; s), (15) 

JO J-oo 

where 

K e (<l>, T, <p; s) = (0, e\K(F, V] e). (16) 



Equation ( 15]) can be regarded as the basis equation for direct sampling 
of CPS distributions from the difference- count statistics in balanced homo- 
dyning. Let us consider the sampling function K e (<ft, J 7 , <p; s) in more detail. 
Expanding the \<j>,e) in photon-number states [see Eq. (||)], Eq. ( |16D can be 
rewritten as 

CO oo 

K e {<f>,F,¥>;s) = (l-e~ 2e )E E K nm {T^-s)e^- m ^e-< n+m \ (17) 

n=0 m=0 



where 



K nm [T, <p; -s) = (n\k{JF,ip\ -s)\m) (18) 



is the sampling function for determining the density matrix in the photon- 
number basis, which has been studied in detail in a number of papers 



H, gg, |T], H, H, 0. It is given by 

K nm {T, V ; -s) = e^ n - m ^f nm (x; s), (19) 

x = J 7 / (^/2\F\), where f nm (x;s) may be expressed in terms of parabolic 
cylinder functions. From Eqs. ( |17D and (|T9| ) we see that K e ((p, J 7 , <p>\ s) is a 
27r-periodic function of the sum phase <p + (p. Since the relation f nm (x; s) = 
fmn(x; s) is valid, K e (<p-> V 9 ! s ) is an even function of the sum phase. For rj 
> 1/2 the function f nm (x, s) is bonded. For large values of x it behaves like 
x -\ n - m \~ 2 , so that K^J 7 ,^; s) is bounded for e > 0. In Figs. |(a) - |(d) 
K e ((f)j F i V 9 ) s ) i s shown for various values of the smoothing parameter e and 
the detection efficiency 77. 

Let us first restrict attention to perfect detection (77 = 1, i.e., s = 0), 
Figs. £|(a) - f|(c). We see that with increasing e the fine structure in the 
dependence on T of K e (<fii 3~ '1 V 9 ! s ) is l° s t- The fine structure obviously results 
from highly oscillating terms in Eq. (p~T|) , i.e., from terms with n,m^> 1. 
These terms become more and more suppressed when e is increased. On the 
other hand, the width of the interval of T in which K t (<j), J 7 , (p; s) is essentially 
nonzero slowly decreases with increasing e. Apart from finenesses, the CPS 
distributions may therefore be expected to reflect typical properties of the 
London phase state distribution even for values of e that are not extremely 
small. Clearly, for e — > the calculation of K e ((p, J 7 , ip; s) becomes an effort, 
because the number of terms in Eq. ( |T7D drastically increases. From Fig.|](d) 
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we see that in the case of non-perfect detection (77 <1, i.e., s<0) the function 
K e ((f), J 7 , if] s) becomes highly structurized, so that increasing effort has to be 
made to sample CPS distributions from the field-strength distributions. 

In general, the CPS distributions p(<p, e) are determined by the field- 
strength distributions p[T ^ <p) in more or less extended intervals of T and 
(p. In particular, Figs. |](a) - (c) clearly show that the intervals of T and 
ip + cf) in which K e (<ft, JF, tp; s) is essentially nonzero are not concentrated 
in the vicinity of the points T = and ip + <fi = tt/2, respectively. Hence, 
the CPS distributions (and the London phase state distribution) cannot be 
obtained from the field-strength distributions at JF = and <y9 + = 7r/2. This 
contradicts the operational phase distribution p(<p) ocp(JF= 0,(p = 7r/2 — 0) 



introduced by Vogel and Schleich [FS 



4 Direct sampling of coherent phase state 
distributions 

As already mentioned, Eq. (|T5| ) can be used for direct sampling of CPS distri- 
butions from the field-strength distributions measured in balanced homodyne 
detection. To give an example, let us consider the determination of CPS dis- 
tributions for a squeezed vacuum state. The result of a computer simulation 
of measurement of a set of field-strength distributions for a squeezed vacuum 
state with mean photon number (n) = 1 is shown in Fig. [5|. In the simula- 
tion perfect detection is considered and 30xl0 4 events at 30 phase values are 
assumed to be recorded. The feasibilty of direct sampling of CPS distribu- 
tions p((f), e) from the measured field-strength distributions is demonstrated 
in Fig. ^| for e = 0.1. 

Comparing the sampled CPS distribution with the calculated distribu- 
tion, we see that the measured double-peak structure is in good agreement 
with the theoretical prediction. In particular, in the intervals of phase in 
which reduced phase fluctuations are observed a sufficiently well agreement 
between theory and experiment is found. Relative large discrepancies be- 
tween the measured and the calculated curves are seen in the intervals of 
phase in which enhanced phase fluctuations are observed. Clearly a bet- 
ter agreement and a higher accuracy can be achieved when the number of 
phases at which the field-strength distributions are measured and the number 
of events recorded at each phase are increased. 



5 Summary and concluding remarks 

We have studied the integral relation of the CPS distributions p((fi, e) to 
the phase-parameterized field-strength p(J- ', (p) measurable in balanced ho- 
modyne detection and calculated the corresponding integral kernel. Since 
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the CPS distributions can be regarded as smoothed London phase state dis- 
tributions, the degree of smoothing being determined by the parameter e in 
the coherent phase states \z = e~ e e % ^) (e>0), in the limit e — > the London 
phase state distribution is observed. 

With decreasing value of the smoothing parameter e the information 
about the phase statistics is increased, because the smoothing parameter 
e controls the interval of the field strength used for "probing" the phase 
statistics of the state under consideration. With decreasing value of e this 
interval is increased and comprises the whole field-strength axis in the limit 
e — > 0. Hence, with decreasing value of e the determination of the integral 
kernel requires inclusion in the calculation of increasing values of the field 
strength. This point must be considered very carefully when the state un- 
der consideration contains large field-strength values that are desired to be 
included in an analysis of the phase statistics. 

In the paper we have calculated the integral kernel using an expansion in 
the photon-number basis. Since in this basis with decreasing value of e the 
values of the photon numbers which must be taken into account increase, 
highly oscillating functions are involved in the expansion and much effort 
must be made to obtain results with reasonable accuracy. This difficulty 
might be overcome by using consequently a field-strength basis and avoiding 
the detour via the photon-number basis, as has recently been demonstrated 



for direct sampling of the density matrix in a field-strength basis [45 



The integral relation between the CPS distributions and the field-strength 
distributions can be used for direct sampling of the CPS distributions from 
the data measured in balanced homodyning, the (bounded) integral kernel 
playing the role of the sampling function. To illustrate the method, we have 
performed a computer simulation of measurement of a set of field-strength 
distributions and sampling from them the CPS distributions for a squeezed 
vacuum state with mean photon number (n) = 1. In this case a smoothing 
parameter of about 0.1 is sufficient to detect the features typical for the phase 
distribution of the state. Using realistic number of phases at which the field- 
strength distribution is measured and realistic numbers of events recorded 
the sampled CPS distributions are found to be in good agreement with the 
calculated ones. 
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Figure 1: The CPS distributions p(<f>, e) for a coherent state with mean photon 
number (n) = 1 are shown for various values of the smoothing parameter e. 
In the limit when e goes to zero p(<f>, e) approaches the London phase state 
distribution (solid line). 
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Figure 2: The CPS distributions p(<f>, e) for a squeezed vacuum state with 
mean photon number (h) = 1 are shown for various values of the smoothing 
parameter e. In the limit when e goes to zero e) approaches the London 
phase state distribution (solid line). 
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Figure 3: The CPS distributions e) for a squeezed vacuum state with 
mean photon number (n) = 2 are shown for various values of the smoothing 
parameter e. Note that compared with Fig. 2 smaller values of e are needed in 
order to obtain CPS distributions that are comparably close to the London 
phase distribution (solid line). 
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(a) (b) 




Figure 4: i^ e (0,jF, ip, ; s) is shown as a function of JF and + for various 
values of the smoothing parameter e and the detection efficiency rj. 

(a) e = 0.1, 7] = 1 (s = 0). 

(b) e = 0.3, 77 = 1 (s = 0). 

(c) e = 0.8, 77 = 1 (s = 0). 

(d) e = 0.1, r7 = 0.8 (s = -0.25). 
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V2\F\p(F^;s) 




Figure 5: The measured set of field-strength distributions p{J- ', ip; s) for a 
squeezed vacuum state with mean photon number (n) = 1 is shown. It is 
obtained from a computer simulation, where perfect detection (i] = 1, i.e., s 
= 0) is considered and 30xl0 4 events at 30 phase values are assumed to be 
recorded (the measured, discrete points are linked by straight lines). 
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Figure 6: The CPS distribution for e = 0.1 which is directly sampled from the 
measured field-strength distributions plotted in Fig. 5 is shown and compared 
with the theoretical curve from Fig. 2 (dashed line). 



17 



